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The Trans-Baikal (or Zabailkal’e) region includes the forest-steppe ecotones south and east of Lake Baikal 
in Russia and has experienced drought for several years. The decline and mortality of birch ( Betula pen- 
dula) stands within the forest-steppe ecotone Trans-Baikal region was studied based on a temporal series 
of satellite data, ground measurements, and tree ring analysis. During the first decade of the 21st century 
birch stands decline and mortality were observed on about 5% of the total area of stands within our 
1250 km 2 study area. Birch forest decline and mortality occurs mainly at the margins of stands, within 
the forest-steppe ecotone on slopes with direct insolation. During the first decade of the 21st century 
summer (June-August) precipitation was about 25% below normal. Soil water content measurements 
were lowest within dead stands and highest within healthy stands and intermediate within damaged 
stands. Drought impact on stands was amplified by an increase in summer air temperatures (+0.9 °C) 
in comparison with the previous decade. Tree ring data of “surviving” and “dead” tree groups showed 
a positive correlation with summer/annual precipitation and negative correlation with summer air tem- 
peratures. Temperature and precipitation extreme anomalies tend to occur in the region with a period of 
about 27 years. The observed anomaly was the most severe since the beginning of meteorological obser- 
vations in the year 1900. Data for the other sites showed a positive climate impact on the growth and 
expansion of Siberian forests. That is, the same species ( B . pendula ) showed considerable increase (1.4 
times both in height and stem volume) during 20th-21st centuries as temperature increased but precip- 
itation remained at adequate levels. 

© 2012 Elsevier B.V. All rights reserved. 


1. Introduction 

Drought-induced forest decline and die-off during the last dec- 
ades has been documented on every continent (Allen et al., 2009). 
Forest mortality caused by elevated temperatures and water stress 
was reported for Europe, including increased mortality among 
many tree species in Spain (Penuelas et al., 2001), France (Breda 
et al., 2006; Landmann et al., 2006), and increases in mortality of 
Pinus sylvestris in Switzerland and Italy (Bigler et al., 2006; Dobber- 
tin and Rigling, 2006), and Populus tremuloides mortality across a 
million hectares was observed in Canada (Hogg et al., 2008) and 
in the US (Anderegg, 2012). Drought-induced tree mortality was 
also documented for over a million hectares of Pinus edulis in the 
south-western part of USA (Swetnam and Betancourt, 1998; Bres- 
hears et al., 2005; Raffa et al., 2008; van Mantgem et al., 2009). In 
Russia, birch stands decline was described for southeast Siberia by 
Anenkhonov (2008). 

Future climate changes are likely to include further increases in 
drying, frequency and severity of droughts in some forested areas 
(Christensen et al., 2007; Seager et al., 2007; Sterl et al., 2008; Allen 
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et al., 2009). Negative impacts on forests include reduced growth 
and increases in stress and mortality caused by synergy of air 
temperature and precipitation change and climate-induced 
changes in the dynamics of dendrophyl insects and fungi (Lucht 
et al., 2006; Scholze et al., 2006; Lloyd and Bunn, 2007). 

Forested areas are especially sensitive to climate variables with- 
in ecotones, the transitions between different vegetation types, 
where tree growth is limited by temperature and/or precipitation 
(Holtmeier, 2009). Climate change scenarios predict aridization 
and risk of conversion of forest-steppe ecotones into steppe or de- 
serts, and northward expansion of steppes (IPCC, 2007). 

The goal of this paper is an analysis of climate-induced impacts 
on birch stands within the forest-steppe ecotone in the Trans- 
Baikal area, the transition between boreal forests and Mongolian 
steppes and deserts (Fig. 1). We hypothesize that air temperatures 
and precipitation changes during past decades will impact tree 
vegetation vigor within the Trans-Baikal forest-steppe. 


2. Study areas 

The study area location is shown in Fig. 1. The Trans-Baikal area 
(also known as Zabailkal’e) is located to the south and east of Lake 
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Fig. 1 . Trans-Baikal study area location (1). 2-5 - show locations for sites of former studies (2 - Western Siberia forest-steppe ecotone, 3 - alpine ecotone in Altay-Sayan 
Mountains; 4 - larch - mixwood taiga ecotone; 5 - northern forest-tundra ecotone; Kharuk et al., 2006-201 1 ). Square indicates on-ground area, Landsat and Quickbird scenes 
location (116°E/ 51°25'N cent point coordinates). 


Baikal in southern Siberia (53°N117°E) and is composed of hills 
and mountains with elevations up to 2500 m. 

Forests were presented by a mosaic of birch-dominated stands 
(Betula pendula ) embedded within the southern Siberian steppe, 
and occupying mainly north facing mountains slopes (up to 
1200 m in elevation). Scots pine ( Pinus silvestris) and Siberian larch 
( Larix sibirica ) were presented by rare specimens. Within study 
area mean density of stands was 1100 trees/ha, mean age was 
about 40 years old; mean tree heights, dbh, and stem volume were 
12 m, 13.5 cm, and 90 m 3 /ha, respectively. South facing slopes 
were covered by grass communities. Soils were rocky podzols. 
Ground water level was at about tens meters. 


2.1. Climate 

The climate is severe continental with harsh cold winters and 
warm summers. During the last four decades mean annual air tem- 
perature and precipitation were -2.7 °C and 330 mm. Mean sum- 
mer (June-August) and winter (December-February) temperature 
and precipitation were +16.5 °C, -24.0 °C, and 232 mm and 
8 mm, respectively. During the last two decades (vs the similar 
previous period) an increase in mean annual temperature 
(+1.1 °C; p < 0.01), a decrease in annual (-37 mm, or >10%) and a 
decrease in summer (-38 mm, or >15%; p<0.05) precipitation 
were observed. 
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3. Materials and methods 

3.1. Satellite data 

A temporal series of satellite data were compiled from Terra/ 
MODIS (http://cimss.ssec.wisc.edu/dbs/China201 1 /Day2 /Lectures/ 
MODIS_MOD13_NDVI_referenc.pdf) and Landsat scenes (https:// 
wist.echo.nasa.gov; http://edcsnsl7.cr.usgs.gov). Ancillary data 
sets included MODIS-based vegetation map with 250 m resolution 
(TerraNorte, http://terranorte.iki.rssi.ru), topography data (http:// 
www.cgiar-csi.org), and a Quickbird scene (~2 m resolution). Ter- 
ra/MODIS and Landsat data were used to generate temporal row of 
vegetation indexes (Enhanced Vegetation Index, EVI, and Normal- 
ized Vegetation Index, NDVI) analysis within the study area. EVI 
was defined as: 

EVI = G ■ (p NIR - Prei ) ■ (p MR + Cl ■ p rei - Cl ■ p blue + L)-' 

where p NIR , pred and p bhie are atmospherically-corrected surface 
reflectance in MODIS bands #1 (620-670 nm), #2 (841-876 nm) 
and #3 (459-479 nm); L is the canopy background adjustment; Cl 
and C2 are the coefficients of the aerosol correction; G - gain factor 
(Huete et al., 1999). EVI, as well as NDVI is responsive to canopy 
structural variations, including leaf area index (LAI), canopy type, 
plant physiognomy, and canopy architecture (Huete et al., 2002). 
Analysis revealed that both Vis, EVI and NDVI performed equally 
well (Wardlow et al., 2010). Technically EVI data are available as 
a ready MODIS product MOD13Q1 from NASA’s site (EOSDIS: 
http://reverb.echo.nasa.gov). In this study EVI values covered the 
period from 2000 to 2010. 

The MODIS vegetation map and topography data were used to 
generate a mask of forest stands. Landsat and Quickbird scenes 
were used for detailed analysis where ground measurements were 
acquired within the study area (Fig. 6). Both Landsat scenes were 
taken within the same vegetative period [i.e., 28 July 2000 (Landsat 
5) and 24 July 2010 (Landsat 7)]. Since the sensors were different, 
the Landsat-5 scene was calibrated against Landsat-7. Geometrical 
correction was performed using “rubber sheeting” (Erdas, 2010). 
Radiometric corrections were made for bands #3 (630-690 nm) 
and #4 (760-900 nm) based on spectral signatures of ground sur- 
faces with relatively stable spectral characteristics (sands, water 
bodies, rock outcrops) (Fig. 2). Within the Landsat scenes identifi- 
cation of dead, damaged and healthy stands was based on ground 
based observations and Quickbird data. 



Fig. 2. Landsat-5 against Landsat-7 regression data for channels 3 and 4. Data were 
calibrated based on the 335 reference points within area of about 34,0000 ha. 


Burned areas were excluded by using a mask which was gener- 
ated based on the MODIS Burned Area Product (MCD45 MODIS: 
https://wist.echo.nasa.gov). Wildfires within the study area were 
observed mainly within the steppe communities. So annual forest 
stand damage was <5% of the total burned area mapped by MODIS. 
Burns were regenerated with grasses and restored to the pre-fire 
vegetation index values during the following year (E. Ponomarev, 
personal communication). Based on this, it was assumed that 
MODIS spectral signatures for a given year were affected by cur- 
rent year fires only. 

3.2. Ground studies 

Ground data were acquired during first 10 days of July 2011. 
Forest inventory and soil measurements were made along tran- 
sects from dead into healthy stands (N= 7; mean length = 70 m; 
elevations range from 460 to 870 m). Measurements were made 
within circular (r = 9.8 m) test plots (N = 21 ) along transects. Mea- 
sured parameters were tree heights, dbh (diameter at breast 
height) and tree status (healthy, declined or dead) (Forest Inven- 
tory Instruction for Russian Forest, 1995). 

Soil moisture measurements were based on samples taken at 
the 25 and 50 cm depth. Weather conditions in the area prior to 
and during measurements were dry with no precipitation observed 
during the sampling periods and only 7-10 mm a week before re- 
corded at the nearest weather station. Soil samples (N = 131 ) were 
obtained from 21 soil profiles. Samples were placed into an alumi- 
num box and immediately weighed. Samples were oven dried in 
the laboratory at 85-90 °C and weighed again. Percent of soil mois- 
ture was calculated based on the weight difference measured pre- 
and post-drying. 

3.3. Normalized hydrothermal index 

Stands water supply was estimated based on Normalized 
Hydrothermal Index (NHI): 

NHI = PnormAnorm, 

where p norm = [(p ( - p min )/(p max - p min ) + 1], t norm 

— [(tj — hnin)/(tmax — fmin) “t" IJj 

Pnorm and t n0 rm are the normalized precipitation and temperature, 
correspondingly scaled to range [1,2]; i is the given year; p* is 
sum of summer monthly precipitation [mm], t 2 is mean summer 
(June-August) temperature [K], p min , p max , t min , and t max are mini- 
mum and maximum summer precipitations and temperatures, cor- 
respondingly during 1901-2009 years. Normalization within range 
[1,2] allowed to equalize an input of precipitation and temperature 
into index values. 

NHI index was similar to Selyaninov hydrothermal index 
(Meshcherskaya and Blazhevich, 1997), which was defined as: 

HI = R/0.m, 

where R is the sum of precipitation during a vegetation season 
when t > +10 °C, t is the sum of the temperatures for the same per- 
iod. Kelvin degrees [K] were used to avoid negative values when 
using Celsius scale. Normalized Hydrothermal Index Anomaly 
(NHIA) were calculated based on: 

NHIA = NHIj - NHImean, 

where NHI 2 is NHI for the given year i and NHI mea n is mean NHI for 
the given period. 


388 


V.I. Kharuk et al./ Forest Ecology and Management 289 (2013) 385-392 





Fig. 3. Air temperature (a), precipitation (b) and hydrothermal index (c) anomaly in the Trans-Baikal zone. 1 - annual, 2 - summer data; 3 and 4 - summer and annual data 
filtered by 13-years moving window. 


3.4. Climate data 

Climate data were obtained from KNMI Climate Explorer 
(http://climexp.knmi.nl). CRU TS3.1 provides climate data aver- 
aged for a cell size 0.5° x 0.5° (~70 x 70 km; Mitchell and Jones, 
2005). 

3.5. Dendrochronology analysis 

On each test plot 3-5 trees were cut at the root neck level and 
disks were obtained for dendrochronology analysis. The surface of 
each sampled disk was sanded or scraped to enhance visibility of 
the tree rings. The widths of tree rings were measured with 
0.01 mm precision using a linear table instrument (LINTAB-III). 
Samples were put into two groups: “dead” (22 specimens) and 
“survived” (11 specimens). The mean tree ages of “dead” and “sur- 
vived” groups were 54 and 56 years, respectively. The TSAP and 
COFECHA computer programs were used in tree ring analysis 
(Holmes, 1983; Rinn, 1996). 

4. Results 

4.1. Climate 

Analysis of available climate data (since 1900) showed positive 
trends during 1960-2009 for annual, summer (June-August), and 
winter air temperatures ( p < 0.05), and for summer temperatures 
during 1928-1947 (p < 0.01; Fig. 3). A negative summer tempera- 
ture trend was observed during 1947-1960. Annual and summer 
precipitations showed positive trends in 1965-1990 (p<0.01), 
and negative trends during 1957-1970 (annual and summer pre- 
cipitation; p < 0.05) and 1988-2006 years (summer precipitations; 
p < 0.01). The trend of summer hydrothermal index was negative 
during 1930-1946, 1957-1970 and in 1983-2009 years 

(p < 0.01), and positive during 1970-1983 (p < 0.05). Fourier-anal- 
ysis revealed a ~27 years cycle in summer air temperature and 
precipitation. The processed data (filtered by 13-years filter, i.e., 
~l/2 of the 27-year detected cycle) are presented on Fig. 3. 

4.2. Trends of satellite data 

Satellite-derived vegetation indexes of birch stands showed a 
decrease during the 2000-2010 time period. For about 5% of the to- 
tal area of stand’s (or about 1250 km 2 ) enhanced vegetation index 
(EVI) decrease was significant at p < 0.05 (Fig. 4). Landsat-derived 



Year 

Fig. 4. MODIS-derived EVI values for the birch stands. Trend were significant at 
p < 0.05. Observed effect of EVI decrease is about two times higher (13% per decade 
vs max. 7% per decade) than MODIS sensor degradation. 



Fig. 5. Birch stand’s NDVI values for the years 2000 (1 ) and 2010 (2); 3 - delta NDVI 
(NDVI 2 oio - NDVI 2 ooo)- Data are significant for dead and damaged stands (p < 0.01). 
Note: all stands used for the study were alive in the year 2000. 

data showed NDVI values decreased with the maximum rate of de- 
crease occurring for dead stands (Fig. 5). 
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Fig. 6. MODIS-derived sketch-map classification: location of the drought damaged and dead stands shown by symbols. Total affected area (~40,0000 km 2 ) bounded by the 
dotted line. Insets. Left: Landsat scene fragment which was used for NDVI-based stands vigor classification (Fig. 5); triangles indicated transect location. Right: Quickbird 
scene showed healthy (1), damaged (2) and dead (3) stands location. 


4.3. Ground data 

Ground based measurements showed that birch forest decline 
and mortality was mainly observed at the margins of stands within 
the forest-steppe contact zone. Drought damaged and dead stands 
occurred in strips of forests that were about 22 ± 5 m and 
40 ± 1 1 m wide, respectively (Fig. 6). Damaged stands were located 
on the sunlit slopes (i.e., south, southwest and southeast facing 
slopes). An example of the observed damage is shown in the inset 
photograph of Fig. 7). Soil water content was found be minimal 


within dead stands, maximum within healthy stands and interme- 
diate within declining stands (Fig. 7). There were no signs of insect 
or fungi attack on the trees. 

Birch stand mortality promotes the rapid establishment of a 
grass cover due to a decrease in competition for water, nutrients 
and light. Height of grasses established within damaged forest area 
was about 2-3 times taller than within drought unaffected stands 
(30-40 cm vs 10-15 cm). 

Dendrochronology data for two tree groups (“dead” and “sur- 
viving”) showed that incremental growth values were similar up 
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Fig. 7. Soil moisture content at 25 and 50 cm depth for the healthy, damaged and 
dead stands. Inset: Birch dieback pattern shown occurs in strips at the margins of 
forest stands. 

to the beginning of the 1970s. Later the stands that eventually died 
(dead group) showed lower incremental growth values (Fig. 8). The 
majority of tree dieback occurred during 2003-2005 (Fig. 8). 

Tree ring widths of “surviving” and “dead” groups were corre- 
lated with climate variables. For “survivors” this relation was de- 
scribed by: 

T w - —0.30 t ma 0.63 P m ai R = 0.70, 

where T w is average tree ring width, t ma and Tma are May-August 
temperature and precipitation, respectively. The equation is signif- 
icant at p< 0.05. Coefficients for t ma and Tma are significant at 
p < 0.01 and at p< 0.001, respectively.For the “dead” group the 
equation is as follows: 


T w = -0.42 t ma + 0.44 P ma ; R 2 = 0.55; equation significant at 
p < 0.03. The coefficients for t ma and P ma are significant at p < 0.01 
and at p < 0.005, respectively. 


5. Discussion 

The Trans-Baikal forest-steppe ecotone is an area of low precip- 
itation (about, 330 mm total). Summer precipitation decreased 
(about ~ 25%) during the first decade of the 21st century (vs similar 
previous period; p < 0.01 ; Fig. 3B). This precipitation decrease coin- 
cided with an observed decrease in vegetation indices (EVI, NDVI) 
for an area of about 1250 km 2 (i.e., area of birch stands where sig- 
nificant decrease of EVI was recorded), or 5% of the total area of 
birch stands; Figs. 4, 5 and 6); both events were correlated 
(r= 0.76). It should be noted that the Terra/MODIS sensor showed 
degradation over the time period of this investigation, i.e. decreas- 
ing of the recorded radiometric response at the rate of about from 
1-2% to 4-7% per decade, depending on the spectral channel 
(Wang et al., 2012). This sensitivity decrease may cause a temporal 
decrease in the EVI values with maximum estimated effect of 
about 7% per decade. Meanwhile, the observed decrease of EVI 
was about twice higher (13% per decade), i.e., the effect of birch 
stands EVI decrease was not an artifact. 

The birch stands studied were strongly sensitive to summer 
precipitation anomalies because of extremely low winter precipi- 
tation (<10 mm), low fall (~50 mm) and spring (~40 mm) precip- 
itation. In addition, soils (rocky podzols) within the study area 
have a low water holding capacity, with water table location far 
below the birch root zone (i.e., within 10-100 m, whereas depth 
of birch roots (at Age = 50 years) is about 4.0 m; Vorobjeva, 
1959). The exception is a minor portion of the region underlain 
by permafrost, where available moisture was higher due to thaw- 
ing permafrost, as well as in areas of perched water tables due to 
permafrost. Soil water content within dead stands was about 30- 
50% lower than within alive stands (Fig. 7). Stands mortality occurs 
in a strip pattern located within the birch-steppe boundary on the 
southern and south-facing slopes with relatively drier and warmer 
conditions. A similar observation was described for aspen stands 
mortality in the western US (Huang and Anderegg, 2012). Drought 



Fig. 8. Dynamics of tree ring width for the “dead” (1) and “survived” (2) stands; 3 - confidence level (p < 0.05), 4 - percentage of sampled dead trees (N = 22) dying in the 
given year. Data were filtered with a 3-years moving window. 
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impact on stands was amplified by an increase in summer air tem- 
perature (+0.9 °C in comparison with the previous decade; 
p<0.02; Fig. 2A). Recent advances in the understanding of tree 
mortality mechanisms suggest that forests could be particularly 
sensitive to increases in temperature in addition to drought alone 
(Allen et al., 2009). The hydrothermal index, which includes both 
temperature and precipitation data, showed the greatest decreases 
have occurred since the year 2000 (Fig. 3c). 

All drought-killed trees died in the period of 2000-2010 with 
maximum mortality occurring during 2002-2004 years, i.e., within 
the maximum period of drought (Fig. 3b and c). Notably, the initial 
tree ring increment for the “dead” and “survived” groups was sim- 
ilar. A deviation occurred between groups at the end of 1970s 
through the beginning of the 1980s, i.e., during the previous 
drought cycle that peaked in 1972 (Fig. 8). With precipitation 
and hydrology index increases during the late 1980s tree ring 
growth also increased. With ending of the current drought 
(2009-2010) “survivors” again showed a growth increase (Fig. 2b 
and c and Fig. 8). The last event was accompanied by leveling off 
of EVI following the decline from 2002-2008 (Fig. 4). Tree-ring 
studies of aspen growth under drought also document large de- 
creases in basal trunk growth during drought years (Hogg et al., 

2005) . 

Along with temperature and precipitation, wildfires may have a 
potential input to the forest decline and dieback. However, wild- 
fires occurred mainly within the steppe grass communities. 
Ground studies showed that burns were infrequent and wildfires 
were not the primary cause of birch mortality. Similar observations 
were reported by Anenkhonov (2008). 

The 27 year cycles in precipitation and air temperature suggest 
that similar birch die-off occurred in the past and may occur again 
in the future. The next drought can be expected around the year 
2030 and expected to be more severe since even conservative cli- 
mate scenarios predict an increase of extreme hot temperatures 
(IPCC, 2007). This is supported by the data obtained. Based on 
the results temperature and hydrology index anomalies associated 
with the birch decline were extreme during period of meteorolog- 
ical observations (Fig. 2a and c). Thus, increasing birch stand dam- 
age may be expected in the future. In the light of this fact larch 
should be considered a preferable species for re-forestation since 
larch exceeds birch in drought resistance. Moreover, according to 
local historical data larch was in the past a dominant species in 
the Trans-Baikal area but was selectively removed by harvesting. 

Finally, the described birch die-off phenomenon in the Trans- 
Baikal region is a part of the broader climate-induced changes 
within Siberian forests. Published data showed a positive climate 
change impact on the growth and expansion of Siberian forests. 
Thus, similar species ( B . pendula) in West Siberia forest-steppe eco- 
tone under higher temperatures, but sufficient precipitation, 
showed considerable increase (1.4 times both in height and stem 
volume) during 20th-21st centuries (e.g., site 2 in Fig. 1). Cli- 
mate-induced forest density increases and upslope tree migration 
were reported for an alpine forest-tundra ecotone in Altai-Sayan 
Mountains (site 3 in Fig. 1, Kharuk et al., 2009, 2010, 2011). Tree 
migration into tundra and stand densification was observed within 
northern forest boundary in Siberia (Fig. 1, site 4; Kharuk et al., 

2006) and in the Polar Ural Mountains (Devi et al., 2008). Along 
with this, climate-driven migration of Siberian pine and spruce 
into traditional larch habitat has been reported (Fig. 1, site 5; Khar- 
uk et al., 2007). 
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